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Abstract 

We study possible contribution of the Majorana neutrino mass eigenstate I'h dominated by a 
sterile neutrino component to neutrinoless double beta (Oz^/3/3) decay. From the current exper- 
imental lower bound on the Oz^/3/3-decay half-life of ''^Ge we derive stringent constraints on the 
i^h — Ve mixing in a wide region of the values of mass. We discuss cosmological and astrophysical 
status of Uh in this mass region. 
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The sterile right-handed neutrinos uji, not participating in the electroweak interactions, 
are natural candidates for the extension of the standard model (SM) field contents. They are 
very typical for various scenarios of the physics beyond the SM. It has long been recognized 
that the sterile neutrinos may have plenty of phenomenological implications. One of the 
celebrated examples is the seesaw generation of tiny Major ana neutrino masses for the three 
observable very light neutrinos via a very large Majorana mass term of the right-handed 
neutrinos. In this case together with the very light neutrinos there also appear very heavy 
Majorana neutrino mass states with the typical masses of ~ 10^^ GeV. However, in a more 
general case of the sterile-active neutrino mixing there may also appear additional neutrino 
states with arbitrary masses m/j. The neutrino mass eigenstates u^, dominated by the 
sterile z//j neutrino weak eigenstates, contain some admixture of the active i'e,fi,T neutrino 
species that allows the Uh to contribute to various processes, in particular, to those which 
are forbidden in the SM by the Lepton Number or Lepton Flavor Violation (LNV or LFV). 
They may also modify the interpretation of cosmological and astrophysical observations. 
Therefore, the masses rrih of Uh neutrino states and their mixing with the active neutrinos 
are subject to various experimental as well as cosmological and astrophysical constraints. 
Various implications of the sterile neutrinos have been extensively studied in the literature 
within different scenarios of neutrino mass generation 

Experimentally, the neutrino states can be searched for as peaks in differential rates 
of various processes and by direct production of Uh followed by their decays in detector (for 
summary see Ref. j^). The z/^ can also give rise to a significant enhancement of the total 
rates of certain processes if their masses, rrih, have the values in appropriate regions specific 
for each of these processes 4 , ^. This effect would be especially pronounced in reactions that 
are forbidden in the SM. The LNV and/or LFV processes belong to this category. Many 
of them are stringently restricted by experiment and allow one to derive stringent limits 
on the Uh contribution. Unprecedented sensitivities have been reached in the experiments 
searching for neutrinoless nuclear double beta decay (Oz//5/5-decay) j^. The current hmits 
on the half-life of Oi//5/?-decay are known to stringently constrain various parameters of the 
LNV physics, in particular, those which correspond to the Majorana neutrino contributions. 

The sterile neutrino contribution to Oz//?/?-decay has been first studied in the pioneer 
work Ref. P| within a specific model with the two sterile neutrinos mixed only with the z/g 
flavor state. In the present paper we develop this analysis by a more careful treatment of 
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the nuclear Oz//5/3-decay matrix elements on the basis of a well motivated nuclear structure 
approach known as the renormalized proton-neutron Quasiparticle Random Phase Approx- 
imation(QRPA) 0]. Also, we adopt a general phenomenological approach to the neutrino 
mass generation without referring to a specific model. We consider an extension of the SM 
with the three left-handed weak doublet neutrinos v'j^^ = {v'lg^-, v'l^-, Vi^) mixed with n species 
of the SM singlet right-handed neutrinos z/^j = {v'jn-, ...i^^^). In general, in this scenario there 
appear + 3 Majorana neutrino mass eigenstates Vi = Ul^v'^^, where Uki is the correspond- 
ing mixing matrix. For simplicity, we assume that in addition to the three conventional 
light neutrinos there exists only one Majorana neutrino mass eigenstate Vh-, dominated by 
the sterile neutrino species, with an arbitrary mass m/j, which may mix with all the active 
neutrino weak eigenstates, l'e,^l,T■ We study possible contribution of this Vh neutrino state 
to Oz//3/5-decay via a nonzero admixture of Ve weak eigenstate. From the non-observation of 
this LNV process we derive stringent limits on the Vh ~ mixing matrix element Veh in a 
wide region of the values of rrih. 

The contribution of Vh neutrino state to the Oi//5/3-decay amplitude is described by the 
standard Majorana neutrino exchange between the two /^-decaying neutrons. The corre- 
sponding half-life formula for a transition to the ground state of final nucleus is given by the 
expression 



Here, mg is the mass of electron and Gqi is the integrated kinematical factor (for its numerical 
values see e.g. Ref. j^]). The nuclear matrix element M^^(rnh) we write in the form 

M°"(m;,) = (/f^(m;„ri2)l + i/Gr(mft,ri2)ai2 + i^T(m;„ri2)Si2) (2) 

with the three terms corresponding to the Fermi, Gamow- Teller, and tensor contribu- 
tions, which depend on the neutrino mass vrih and on the relative coordinates defined as 
ri2 = I"! — i'2, f\2 = |ri2|, ^12 = i"i2/ti2, whcrc ri and r2 are the coordinates of beta decay- 
ing nucleons. We also introduced the operators: 5*12 = 3((Ti-ri2)(cT2-ri2) — o"i2j o"i2 = f?i-'?2- 
The neutrino-exchange potentials HK{iTLh,ri2), with K = F,GT,T, have the following 
form: 
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where fxigru) = joiqru), joigr^) and j2{qri2) for K = F, GT and T, respectively (jo and 
j2 are spherical bessel functions). E^ ^ , ^ , and EJ- are the energies of initial, final, and 
intermediate nuclear states, respectively. The mean nuclear radius is given hy R = tqA^^^ 
with Tq = 1.1 fm, and = 1.25 is the axial nucleon constant. The functions hxiq"^) are 

hriq ) - -9y{q ), hr{q )-- 3 4m2 + 3 Am^ ' 

hcAq ) = gAq ) - 3 + 3 4^2 + 3 Am^ ' 

For the nucleon form factors we use the standard parameterizations: 

gviq = (1 + g VA^)-', gAiq ') = ^?a(1 + q V ^lY^ dniq = (/ip - /in)(7y(g 

^p(g 2) = 2mp^^(g ^){q ^ + ml)-\ where (/ip-/i„) = 3.70, Ay = 0.84 GeV, = 1.09 GeV, 

and m^r is the pion mass. 

We calculated the nuclear matrix element in Eq. Q for the transition 
™Ge(0^^ ) -^^^Se(0^g ) within the proton-neutron renormalized QRPA 0. We analyzed 
the single-particle model space which involves the full 2 — Ahu major shells both for pro- 
tons and neutrons. The single particle energies were obtained using the Coulomb-corrected 
Woods-Saxon potential. We calculated the two-body G-matrix elements from the Bonn 
one-boson exchange potential on the basis of the Brueckner theory. Since the considered 
model space is finite the pairing interactions have been adjusted to fit the empirical pairing 
gaps. In addition, we renormalized the particle-particle and particle-hole channels of the 



G-matrix interaction of the nuclear Hamiltonian following the procedure outlined in Ref. 
In this way we obtained the nuclear matrix element with a weak dependence on the nuclear 
structure parameter space. Finally, in order to take into account the two-nucleon correla- 
tions we multiplied the transition operators by the square of the correlation Jastrow-like 
function 

Having the nuclear matrix element M^'^{mh) calculated, we can derive the Oz//3/3 decay 
limits on the mass rrih of the neutrino and its mixing Ueh with the neutrino weak 
eigenstate. Applying the presently best lower bound on the Oz//3/3-decay half-hfe of ^^Ge 

> ^1%"'"" = 1-9 X lO'Vears (5) 

we derive from Eq. the |f/ehP ~ n^'h exclusion plot shown in Fig. ^ In Fig. ^we also 
show typical domains excluded by some other experiments summarized in Ref. j3|. However, 
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these domains are just indicative, because most of the previous bounds were obtained for 
some fixed values of rrih. For convenience, we interpolated this set of experimental points 
by continuous curves in different intervals of m/j. Let us also point out the following. The 
constraints listed in Ref. ^ are based on the searches for peaks in differential rates of various 
processes and the direct production of states followed by their decays in a detector. In 
Ref. y it was argued that in this case the results of data analysis depend on the 
total decay width, including the neutral current decay channels. The latter have not been 
properly taken into account in the derivation of the mentioned experimental constraints. 
However, the neutral current z/^ decay channels introduce the dependence of the final results 
on all of the mixing matrix elements Ueh, U^h and Urh- In this situation one cannot extract 
individual limits for these matrix elements without some additional assumptions, introducing 
a significant uncertainty. In contrast, our Oz//5/5-decay limits involve only Ueh mixing matrix 
element and, therefore, are free of the mentioned uncertainty. This is because in Oz//?/5-decay 
intermediate Majorana neutrinos are always off-mass-shell states and their decay widths are 
irrelevant. 

It is well known that massive neutrinos may have important cosmological and astrophysi- 
cal implications. They are expected to contribute to the mass density of the universe, partic- 
ipate in cosmic structure formation, big-bang nucleosynthesis, supernova explosions, imprint 
themselves in the cosmic microwave background etc. (see, for instance, Refs. 0,|3,Q 
and references therein). This implies certain constrains on the neutrino masses and mix- 
ing, which, however, for unstable neutrinos, such as Uh studied in the present paper, are 
quite weak and involve various uncertainties. The Big-bang nucleosynthesis(BBN) and the 
SN 1987A neutrino signal provide the most stringent and certain constraints [12 [for this 
case. The corresponding constraints on Uh parameters have been derived in Ref. [l^ for the 
restricted mass regions: lOMeV < rrih < 200 MeV for the BBN and lOMeV < rrih < 100 
MeV for the SN 1987A cases. We present the corresponding exclusion plots in Fig. Q for 
comparison with our limits. 

As seen from Fig. ^ our Oz//3/3-decay limits cover much wider region of Uh masses rrih 
and exclude the parts of the |?7e/iP — n^-h parameter space previously unconstrained by the 
laboratory experiments, and the cosmological and astrophysical observations. 

In summary, we examined the possible contribution of the sterile neutrino dominated 
Majorana mass eigenstate Uh to the OujSjS decay. From the currently most stringent lower 
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bound on the 0u(3(3-decaiy half-life of Ge we derived upper limits on the Uh — Ve neutrino 
mixing matrix element |C/efeP in a wide region of the values of m/j, excluding previously 
unconstrained domains in the \Ueh\^ —'rrih plane. Our constraints are free of the uncertainties, 
related to the assumptions on the values of the other neutrino mixing matrix elements U^h 
and Ut^, typical for the other laboratory constraints of this type. 
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Figures 
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FIG. 1: Exclusion plots in the |C/e/iP — "^/i plane. The shaded regions are excluded by the Oz^/?/?- 
dccay [present result], and by the big bang nucleosynthesis and the SN 1987A neutrino observations. 
The dashed line mark typical upper limits from the laboratory searches for massive neutrinos. 
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